Zipping and unzipping of nanoscale carbon structures 
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We demonstrate theoretically that hydrogenation and annealing applied to nanoscale carbon struc- 
tures play a crucial role in determining the final shape of the system. In particular, graphene flakes 
characterized by the linear and non-hydrogenated zigzag or armchair edges have high propensity 
to merge into a bigger flake or a nanotube (the formation of a single carbon-carbon bond lowers 
the total energy of the system by up to 6.22 eV). Conversely, the line of the sp 2 carbon bonds 
(common for pure carbon structures such as graphene or a carbon nanotube) converted into the sp 3 
type by hydrogenation shows an ability to disassemble the original structure by cutting it along the 
line of the modified bonds. These structural transformations provide us with an understanding of 
the behavior of mobile carbon structures in solution and a distinct scenario of how to preserve the 
original structure which would be a crucial issue for their application in carbon-based electronics. 
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INTRODUCTION 

The unique electronic properties of graphene [HQ, in 
particular the high mobility of its charge carriers make 
this recently discovered material an excellent candidate 
for semiconductor electronics [3|, |4[ . The band gap of the 
finite size graphene is defined by the tt states induced by 
the unsaturated dangling bonds at the edges, and there- 
fore, the electronic properties of the finite size graphene 
structures strictly depend on the size and geometry of 
the flakes and also on their atomic edge structure. Since 
these 7r states are spin-polarized, the spin polarization 
of graphene will depend on the atomic edge structure as 
well. It was revealed that under certain conditions (for 
example in an applied electric field [5J or doping of the 
zigzag edges [6]) the spin- up and spin-down states can be 
separated between two zigzag edges such that the half- 
metallicity is achieved. Additionally, if an imbalance in 
concentration of the spin-up and spin-down states are to 
be introduced through edge modification, it would re- 
sult in non-zero magnetization Q. It is then expected 
that the size of the graphene flakes and the atomic struc- 
ture of the zigzag edges are the most crucial parameters 
for application of graphene in semiconductor electronics 
and spintronics. It is already known that the standard 
chemical methods used to fabricate the graphene flakes 
lead to the rough edges [8-10]. Therefore, for the rapidly 
developing field of carbon electronics, the technology of 
"cutting" graphene into nanoscale flakes of controllable 
shape with linear edges is immensely important. Inter- 
estingly, there is also a growing demand for generation 
of large graphene sheets which have already found their 
applications in touch screens and solar cells [Til ] . 

Development of the technology for cutting the carbon 
structures has been under consideration since the last 
decade and began with the discovery of carbon nanotubes 
[l2|. Originally it was observed that cutting and un- 
zipping of nanotubes can be performed in highly acidic 
environment where hydrogen ions supposedly break the 
carbon-carbon bonds, while further annealing restores 
some of the tubular structures [l2|| and restoration would 



be more sufficient at higher temperatures [13|. It was 
also observed that annealing of nanotubes can lead to 
coalescence of the nearest located tubes into the one 
of a larger diameter [14]. Since then, several technolo- 
gies have been developed to unzip nanotubes [HI but 
the most popular ones are chemical methods which use 
oxidizing agents to break the carbon-carbon bonds [l6| . 
and catalytic hydrogenation using the metal nanoparti- 



cles as scissors [17, 



18j. Therefore, the original idea [12| 
that nanotubes can be disintegrated or integrated due to 
the interaction with the atomic hydrogen lost its appeal 
with the discovery of these new chemical methods. Of 
late, it has received a fresh scrutiny [19| after the discov- 
ery of graphene. These authors have used hydrogenation 
of graphene edges to control the nanoribbon width and 
the edge quality within the plasma etching procedure. 

Therefore, in this work we wish to refocus the atten- 
tion on the method of manipulating the carbon structure 
by their interaction with atomic hydrogen, by providing 
insights on the phenomena of merging and cutting pro- 
cesses in carbon structures with the help of the atomic 
hydrogen. We have simulated the transformation (merg- 
ing and unzipping) of the carbon system such as graphene 
flakes and carbon nanotubes in vacuum with the help 
of quantum chemistry methods. The calculations have 
been performed using the spin-polarized density func- 
tional theory with the hybrid functional UB3LYP/6-31G 
available in the Jaguar 7.5 program [20]. For interac- 
tion and bonding of the carbon atoms with hydrogen 
atoms the B3LYP exchange-correlation functional has 
been found to offer sufficient precision level for the geom- 
etry parameters and energetics [21] thus providing a good 
agreement with the results obtained via a more accurate 
level such as the ab initio M0ller-Plesset second-order 
(MP2) method. The 6-31G basis set was found to be 
adequate to describe the alteration of the total energies 
and lattice parameters in the carbon systems. The initial 
structures of the graphene flakes have been generated us- 
ing the standard parameters of graphene honeycomb lat- 
tice (the carbon-carbon bond length was 1.42 A). Dur- 
ing the geometry optimization (UB3LYP/6-31G) the sys- 
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AE=-8.97eV AE = -37.5eV 
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AE = -10.67 eV AE = -33.98 eV 



FIG. 1: Integration of two isolated flakes characterized by the 
sharp and linear edges (top panel presents the case when the 
flakes are stacked by the zigzag edges, while the bottom panel 
is for flakes stacked by the armchair edges) into one flake of 
larger size: a) and d) two initially separated flakes; b) and 
e) formation of the first carbon-carbon bond; c) and f) the 
merged systems. 

terns were allowed to relax in all directions thereby form- 
ing the shape characterized by the lower total energy. For 
all optimization procedures, the transition of the struc- 
ture from an initial state to the optimized geometry (with 
corresponding structural changes of the graphene lattice) 
during the relaxation procedure occurs without any en- 
ergetical barrier. To perform the hydrogenation of the 
graphene flakes the hydrogen atoms have been attached 
to the required carbon atoms and their positions have 
been further optimized. 

ZIPPING PROCESS 

Recent developments in the technology of fabrication 
of carbon systems have demonstrated that the sharpness 
of zigzag edges in graphene can be achieved by Joule 
heating [22j, i.e., by evaporation of carbon atoms whose 
presence at the edges break the linearity of these edges. 
The same effect is also achieved by H2 -plasma etching 
[l9| . These methods result in sharp and highly crys- 
talline edges that, as we mentioned above, would promote 
the application of graphene in semiconductor electronics. 
In this respect we investigate the properties of graphene 
flakes of rectangular shape and found that when their 
edges are linear and non-hydrogenated the interaction 
between the two flakes is so high that several flakes will 
have the tendency to merge into one big structure. This 
result is extremely important and requires considerable 
attention because it might have valuable impact on the 
technology of fabrication of the touch screens and solar 
cells [11] as a way to create a graphene sheet of desirable 
size through flake integration. 

For our investigations we took two graphene flakes with 
linear zigzag or armchair edges and placed them side by 



side, as shown in Fig. [T] (a) and (d). After relaxation 
of the system, two flakes were merged along the stacked 
edges, thereby generating one flake of bigger size. The 
flakes with the linearly cut zigzag edges zip up along 
those zigzag edges (see (a-c) in Fig. [T]) while in the other 
case the zipping occurred between the stacked armchair 
edges (see (d-f) in Fig. [1]). We presented in Fig. [T] three 
major steps of the merging process: the separated flakes, 
origin of the zipping process, and the merged system. 
Zipping begins by moving the flakes close to each other 
by one seam with the formation of the first carbon-carbon 
bond between the two flakes (see (b,e) in Fig. [1]) and af- 
ter that the rest of the system zip up (see (c,f) in Fig. [I]). 
To show that zipping occurs because it is the energeti- 
cally favored process, we calculated the energy AE (see 
its value in Fig. [T]) as the difference between the total 
energy of the original systems (Fig. [T] (a) or (d) for the 
zigzag and armchair edges, respectively) and the total 
energy of the structure taken in the current conforma- 
tion. The energy \AE\ grows as the flakes move close to 
each other (Fig. [T] (b) or (e) for the zigzag and armchair 
edges, respectively) due to the interaction between them, 
while subsequent formation of each carbon-carbon bond 
increases that energy even further. 

Since formation of a large flake through zipping of two 
flakes of smaller size significantly lowers the total energy 
of the system, it implies that this process should occur 
with any flakes having the following characteristics: non- 
hydrogenated linear edges and ability to migrate. The 
evolution of the total energy with zipping of the flakes is 
defined not only by the zipping mechanism itself but by 
other parameters such as the original flake conformation 
and the distance between the flakes. However, to gain in- 
sight on the merging phenomena we need to estimate the 
energy released as each carbon-carbon bond is made. For 
this purpose we calculated the formation energy of the 
flake integration as AE Z = E^ merged ^ — (E\ +^2), where 
E(merged) is the total energy of the system after integra- 
tion, while E\,E2 are the total energies of the initially 
separated first and second flakes, respectively. We found 
that in the case when zipping involves the zigzag edges 
(see (a-c) in Fig. [I]), the formation energy AE Z was -43.99 
eV, while for zipping along the armchair edges (see (d-f) 
in Fig. [T]) AE Z =-51.20 eV. However, because the con- 
formation of the graphene structure is different for the 
initial and final steps, the formation energy AE Z con- 
tains two major components and one of which is the en- 
ergy released with the generation of each carbon-carbon 
bond while the second one is the reorganization energy 
of the lattice, i.e., the energy required to reorganize the 
graphene lattice with a transition of the system from the 
initial to final states. Therefore, to estimate the energy 
released by formation of the single carbon-carbon bond, 
we need to substitute the reorganization energy from the 
total formation energy. We calculated the reorganization 
energy to be -6.63 eV for merging of the flakes through 
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A E=-15.68 eV A E=-24.29 eV A E= -48.76 eV 



FIG. 2: Formation of a carbon nanotube from two graphene 
flakes placed on top of each other with linearly cut and non- 
hydrogenated zigzag edges: a) two separated flakes, b) origin 
of the zipping process, c) formation of the carbon nanotube. 

zipping of the zigzag edges and -5.78 eV for flakes merg- 
ing through zipping of the armchair edges. As a result we 
found that the formation energy of a single bond made 
between the zigzag edges is ~ -6.22 eV. This value is close 
to that of the dissociation energy of the double carbon- 
carbon bonds which is ~6.4 eV [23j. For the armchair 
edges we obtained the formation energy of a single bond 
to be ~ -5.92 eV, but this value can not be very accurate 
because in the case of zipping of the armchair edges the 
situation is more complicated as the armchair edges can 
not zip up completely without introducing some defected 
sites at the end of the line along which two flakes join (see 

Fig.HKf)) 

As it is impossible to merge the armchair edges with- 
out introducing the defected sites and because the zigzag 
edges seems to have a higher tendency to zip up than 
the armchair edges, the zipping along the zigzag edges 
should receive careful consideration for fabrication of 
large graphene films required for touch screens and solar 
cells technology For these applications, the effect 

of the flake size (finite size effects) should be an impor- 
tant concern. Because the lattice reorganization within 
the zipping process occurs mostly at the zipping edges 
the contribution of the reorganization energy to AE Z 
will gradually decrease with increasing length of the flake 
along the zipping edges. A similar conclusion was reached 
for the second energy component (linked with the forma- 
tion of the carbon-carbon bonds) contributing to AE Z , 
i.e. AE Z will be reduced with an increase of the number 
of carbon-carbon bonds that are to be formed. There- 
fore, an increase in flake size, namely the flake length 
along the axis of zipping, would favor an increase of en- 
ergy AE Z released in the process. Growth of the flake 
size in the other direction will not drastically modify the 
energetics of the zipping process. Moreover, the large for- 
mation energy of the carbon-carbon bond (-6.22 eV for 
the zigzag edges) suggests that the flake zipping would 
occur even at non-zero temperatures as this formation 
energy is much larger than the thermal energy (fc^T). 
The last conclusion is in fact, supported by the experi- 
mental observation of zipping of the damaged nanotubes 

30. 

Because generation of each carbon-carbon bond re- 
leases so much energy, placing the flakes with linearly 
cut zigzag edges on top of each other may lead to oc- 



currence of the zipping process at both sides limited by 
the zigzag edges thereby creating a nanotube. We per- 
formed the relaxation of the system presented in Fig. [2] 
(a) and achieved the rolling of the flakes into a nanotube 
(see Fig.[2](d)). At first, the graphene flakes slowly bend 
into a tubular form (see Fig. [2] (c)) and with that bend- 
ing the zigzag edges start to zip up (see Fig. [2] (d)). 
We calculated the energy of nanotube formation to be 
AE Z =-47.86 eV. Because there are twice the carbon- 
carbon bonds to be made to generate a nanotube from 
two flakes, it was expected that the energy released due 
to the tube formation should be much larger than that 
for integration of two flakes into one of bigger size. To 
confirm that we took the same flakes as in Fig. [2] (a) 
and placed them side by side, zipped them along their 
zigzag edges and calculated the formation energy to be 
AE Z =-36.58 eV. As expected, we found that the genera- 
tion of a nanotube is much more favorable energetically in 
comparison to flakes merging into a flake of larger size. 
For flakes presented in Fig. [2] (a) the energy difference 
between the two final shapes (a nanotube in compari- 
son to a flake) achieves the magnitude of -13.86 eV and 
this difference will grow with elongation of the nanotube 
structure. The reorganization energy in the case of the 
nanotube formation is quite small (-1.21 eV) which gives 
us the formation energy of the carbon-carbon bond to be 
-5.83 eV. In contrast, for the flakes merging into a flake of 
larger size, the reorganization energy remains significant 
(-10.66 eV) giving a bond formation energy of -6.48 eV. 

Therefore, if there are two flakes with non- 
hydrogenated linear zigzag edges, depending on certain 
conditions they can form either a flake of larger size or a 
carbon nanotube. Obviously, the main condition control- 
ling the final shape would be the mobility of the flakes 
and their original position relative to each other. There- 
fore, in a solution where the carbon systems presumably 
have higher mobility, the energetically favored shape of 
the carbon structure after annealing is expected to be a 
carbon nanotube. That is the reason for zipping of par- 
tially unzipped carbon nanotubes back to the original 
nanotube when its edges at the cutting front possess the 
dangling bonds [13| . As formation of each carbon-carbon 
bond drastically lowers the total energy of the system, the 
process involving more bond formation would be favored 
energetically. Moreover, in a geometry of the carbon sys- 
tems with less number of carbon atoms at the edges the 
total energy is lowered. In fact, this observation offers an 
understanding of the experimental data on coalescence of 
the carbon nanotubes in solution [14]. For example, we 
estimated that coalescence of two nanotubes (placed one 
behind other thereby forming a longer tube) of the same 
diameter equal to that in Fig. [2] (d) lowers the total en- 
ergy by 51.34 eV. 

Another important understanding is that the absence 
of hydrogenation at the edges makes the graphene flakes 
chemically attractive to each other. If interaction of the 



graphene flake with the surrounding environment can be 
considered as weak [25j (because interaction with most of 
the adsorbates occurs via the van der Waals forces (26|), 
the interaction between the graphene flakes itself is quite 
significant and may lead to their integration. Therefore, 
if we are to prevent the integration of the carbon struc- 
tures for their application in semiconductor electronics, 
they should be immobilized or being hydrogenated along 
the edges. 

UNZIPPING PROCESS 



4 




FIG. 3: (a) Breaking of the carbon-carbon bond as three pairs 
of hydrogen atoms attached to the carbon atoms line up; (b) 
distribution of the interatomic forces along the z and y axes 
of the graphene plane. 



For graphene hydrogenation, the edges are the ones to 
be first hydrogenated because the energy required to at- 
tach a single hydrogen atom to the edges is about 3.0 eV 
less than that for hydrogen attachment to the graphene 
surface. Hydrogenation of graphene is a thermally acti- 
vated process and is rarely found to occur at room tem- 
perature [27H29|. Therefore, it is possible to hydrogenate 
the graphene edges selectively from the rest of the struc- 
ture by controlling the temperature. For hydrogenation 
of the graphene surface, it appears that the conforma- 
tion characterized by the lower total energy would be 
when the nearest-neighbor carbon atoms belonging to 
different sublattices (there are A and B sublattices in 
graphene Q ) are hydrogenated from different sides of the 
graphene plane [30|. In this case the chemisorption en- 
ergy is characterized by the lowest value and if the surface 
of graphene is fully hydrogenated, the system will turn 
into graphane in its chair conformation 31, 32| . 

For partial hydrogenation of the graphene surface, if 
the first hydrogen atom is randomly attached to the 
graphene surface, the second atomic hydrogen will be 
attached to the nearest-neighbor carbon atom as more 
energetically favorable position [33|. Let us consider the 
process leading to the lowest state, i.e., to graphane in 
the chair conformation that occurs when hydrogenation 
is applied to both sides of the graphene plane. In this 
case, attaching the second hydrogen atom to the nearest- 
neighbor carbon atom with higher probability would oc- 
cur at other side of the plane [30]. To continue the hy- 
drogenation process we have to take into account the fact 
that the carbon atoms located closer to the already hy- 
drogenated part requires less energy to be hydrogenated 
as this energy increases with increasing distance from 
the already hydrogenated carbon atoms [33|. Because 
the pair of hydrogen atoms attached to the nearest- 
neighbor carbon atoms weakens the neighboring bonds, 
under some conditions of hydrogenation the attachment 
of hydrogen atoms may occur linearly (when minor dif- 
ference in magnitude of the binding energy tends to con- 
trol the hydrogenation process). As hydrogenation is a 
thermally activated process 27l-[29|. such a control can 
be obtained at a lower temperature, while its substantial 
increase would actually lead to full hydrogenation of the 




FIG. 4: a) Unzipping of a flake when each carbon atom 
at the edges (including the C1-4 atoms) is terminated by a 
single hydrogen atom, b) Unzipping of a flake for the case 
when C1-4 carbon atoms are terminated by two hydrogen 
atoms before unzipping. 



graphene surface, i.e., to graphane-like structures. This 
effect of weakening of the neighboring bonds by an expo- 
sure to active adsorbates is used for longitudinal cutting 
of the carbon nanotubes by the oxidative process 16|, l34j 
and has been predicted to be the case in hydrogenation 
of the nanotubes [HI . 

For three or more pairs of the hydrogenated carbon 
bonds forming a line as presented in Fig. [3] we found an 
interesting effect: The middle carbon-carbon bond from 
the line of hydrogenated atoms is being broken. To un- 
derstand this behavior, we have to analyze the changes 
brought into the graphene lattice by hydrogenation. If 
two carbon atoms connected by the carbon-carbon bond 
are hydrogenated, it alters the hybridization of this bond 
from sp 2 to sp 3 , thereby enlarging its length from 1.42 A 
to 1.52 A [30]. As a result, the strength of such a carbon- 
carbon bond is significantly lowered. Moreover, hydro- 
genation causes the lattice distortion such that the car- 
bon atoms hydrogenated from different sides of graphene 
are shifted out of its plane in opposite directions (in our 
present work it occurs along the z - axis in Fig. E9 (b)). 
In the ideal honeycomb lattice the inter-atomic forces be- 
tween the carbon atoms are almost compensated within 
the graphene plane (in this work the plane is formed in 
the x and y axes) and only the carbon atoms at the edges 
experience some compressive stress [36[. However, the 
shift of the hydrogenated carbon atoms out of the plane 
leads to the appearance of the uncompensated inter- 
atomic forces arising on the modified bond and those 
forces are pointed out of the plane. The forces function 
as the tensile stress applied to the sp 3 hybridized carbon- 
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carbon bond (along the z and y axes) thus additionally 
reducing its strength. When there are three pairs of hy- 
drogen atoms attached to the carbon atoms as presented 
in Fig. [3] (a), the tensile stress arising at the modified car- 
bon bonds is actually weaker for the first and last pairs 
due to the compensation effect induced by the unmod- 
ified lattice around it. For the middle pair, the tensile 
stress is able to break the carbon-carbon bond. 

We calculated the inter-atomic forces for the ideal 
graphene lattice and found that they are zero in the 
z-direction, while in the x and y-axes they are almost 
compensated within the graphene plane (the remain- 
ing inter-atomic forces is ~ ± 0.009 eV/Bohr of mag- 
nitude). Hydrogenation performed along the line as pre- 
sented in Fig. [3] (a) in the ^/-direction increases the ex- 
isting forces on the hydrogenated carbon atoms up to 
~ ± 0.5 eV/Bohr, while in the z-direction it causes an 
appearance of the inter-atomic forces of magnitude ~ ± 
0.9 eV/Bohr (shown as F y and F z in Fig. 02(b)). There- 
fore, these forces in the z and ?/-axes are larger almost by 
two orders of magnitude than that in the ideal graphene 
lattice and the resulting tensile strain applied to the sp 3 
hybridized bond (see F in Fig. [21(b)) is ~ ± 1.0 eV/Bohr 
from each side. Such a tensile strain appears to be large 
enough to break the weak sp 3 hybridized carbon-carbon 
bond. 

If we continue the line of hydrogenation presented in 
Fig. [3] (a) to the edges, will the tensile stress disinte- 
grate the flake? We observed that the flake separation 
occurs only in the middle of the hydrogenation line (see 
Fig. [4](a)) but not at the armchair edges (see the Ci_4 
carbon atoms in Fig. [4]). Away from the edges the car- 
bon atoms are connected to three nearest-neighbors and 
therefore, after their bonding with the hydrogen atom 
all their valence electrons are saturated thus converting 
the carbon-carbon bond into the sp 3 hybridization type. 
However, the carbon atoms at the armchair edges have 
only two neighbors and their bonding with one hydro- 
gen atom turns the C1-C2 and C3-C4 bonds into sp 2 
hybridization type (the bond length is 1.43 A) instead 
of the required sp 3 type. The termination of the C1-4 
carbon atoms by the two hydrogen atoms (see Fig. 21(b)) 
does not break the C1-C2 and C3-C4 bonds yet but lead 
only to their elongation (the bond length is 1.58 A). To 
break these bonds the termination of the Ci_4 carbon 
atoms (dangling after bond separation) by three hydro- 
gen atoms is required which can not be done within the 
graphene lattice preceding the bond breaking (within the 
graphene lattice those carbon atoms have only two va- 
lence electrons left). 

However, if the line of hydrogenation is shifted as in 
Fig. [5] (a), all the carbon atoms in the line are involved 
into bonding with three nearest-neighbors and attach- 
ment of the one extra hydrogen atom would indeed turn 
all the bonds into the sp 3 hybridization type thereby 
causing their weakening and elongation. The tensile 



strain arising as a result of distortion of the graphene 
lattice breaks all the sp 3 bonds that would lead to dis- 
integration of the original flake (see Fig. [5] (b-c)). We 
calculated the evolution of the total energy (AE) of the 
system as one flake is cut into two flakes. We obtained 
that its disintegration significantly lowers the total en- 
ergy. After the cutting is done (see Fig. [5] (c)), the gen- 
erated flakes take the position supporting the interaction 
of their dangling bonds left on the carbon atoms at the 
zigzag edges. 

Technologically, it is easier to perform hydrogenation 
of the graphene flake from one side of the plane. Tak- 
ing into account the great advantage offered by graphene 
cutting through its hydrogenation, we also investigated 
the behavior of the graphene system in the case of one- 
sided hydrogenation (see Fig.[6|). We took the same flake 
as above and attached the hydrogen atoms at the same 
positions as in Fig. [5] but all from the same side of the 
plane. In the case of one-sided hydrogenation, the hy- 
drogenated carbon atoms are shifted out of the plane 
but in the same direction (see Fig. El (a)). Similar to the 
double-sided hydrogenation, one-sided hydrogenation al- 
ters the hybridization of the carbon-carbon bonds from 
the sp 2 to sp 3 type. However, because the originated 
inter-atomic forces are all pointed in the same direction, 
they do not generate the tensile strain to the bonds but 
initiate the bending of the flake along the line of sp 3 hy- 
bridized carbon-carbon bonds (see Fig. [6] (b)). At some 
point the stress arising at the bend is large enough to 
break these sp 3 carbon-carbon bonds (see Fig. [6] (b)). 
After the bonds are broken, similarly to previous double- 
sided case presented in Fig. 0(c), the disintegrated flakes 
migrate into the position which maximizes the interac- 
tion between the dangling bonds left on the carbon atoms 
at the zigzag edges that lowers the total energy of the 
system. 

Therefore, for the planar graphene flakes if the carbon- 
carbon bonds are turned into the sp 3 hybridization type 
by hydrogenation, it significantly lowers the strength of 
these bonds. Because hydrogenation induces lattice dis- 
tortion thereby inhibiting the planar ity of the flake, it 
causes the occurrence of the stress which is large enough 
to break the sp 3 bonds if they form a line. If two sides of 
the graphene planes are involved in hydrogenation, the 
resultant stress works as the tensile strain applied to the 
sp 3 bonds, while for one-sided hydrogenation, this stress 
initiates the flake bending that eventually breaks the 
bonds. The proposed in-line hydrogenation of graphene 
can be successfully applied to cut the graphene flakes in 
two or more pieces at the same time. Moreover, the same 
procedure can be used to unzip a carbon nanotube into 
a flake characterized by the ideal zigza g ed ges that is in 
agreement with previous investigations [351 ]. 



AE =-8.71 eV 



AE =-19.7 eV 



FIG. 5: The cutting of graphene flakes by formation of a line of carbon atoms being hydrogenated from two sides of the 
graphene plane, a) Shift of the carbon atoms belonging to different sublattices of graphene in opposite direction out of plane 
and elongation of the carbon bond up to ~1.75 A as a result of changing the hybridization of this bond from sp 2 to sp 3 type, 
b) Breaking of carbon-carbon bonds between hydrogenated carbon atoms as a result of bond weakening and bond elongation 
turned on by the sp 3 hybridization and origin of the tensile strain on this bond. Shift of the flakes relative to each other, c) 
Migration of the flakes after their disintegration. 




FIG. 6: The cutting of graphene flakes by formation of a line of carbon-carbon bonds for which atoms are hydrogenated from 
one side of the graphene plane, a) Bending of the lattice due to bonding of carbon atoms with the hydrogen atoms: shift of all 
carbon atom positions in the graphene lattice occurs in the same direction, and elongation of the carbon bond to ~1.75 A as 
a result of changing of the carbon-carbon bond hybridization from sp 2 to sp 3 type, b) Breaking of the carbon-carbon bonds 
between hydrogenated carbon atoms due to the stress induced by modification of the lattice, c) Migration of the flakes after 
their disassembling. 



CONCLUSION 

We have shown that hydrogenation not only plays 
an important role in defining the electronic properties 
of graphene [2], it might also be used to control the 
graphene geometry. We found that flakes characterized 
by the linear and non-hydrogenated edges have the abil- 
ity to merge into a bigger flake or a nanotube through 
zipping of their zigzag or armchair edges. Based on our 
results we suggest that a nanotube would be particularly 
characterized by the lower total energy because it has 
an energetically favorable shape. The energetic prefer- 
ence of a tubular shape comes from the large number of 
carbon-carbon bonds to be made within the zipping pro- 
cess than that for planar geometries. Formation of each 
carbon-carbon bond lowers the total energy of the sys- 
tem by ~ 6 eV. This means that in case of high mobility 
of the carbon structures, to maintain the shape of others 
geometries such as graphene flakes, hydrogenation has to 
be applied to all the edges. However, hydrogenation of 
the graphene flakes has another tricky issue. If its appli- 
cation on the graphene surface creates the line of carbon- 
carbon bonds whose both atoms are hydrogenated, it can 
lead to a separation of the flakes into several pieces. The 



reason for such a behavior is that hydrogenation causes 
the alteration of the bond hybridization from sp 2 to sp 3 
type that significantly suppresses the bond strength and 
the resulting stress turned on by lattice distortion under 
some conditions can be large enough to break those con- 
verted sp 3 bonds. Such cutting of the graphene flakes 
can be performed through one-sided or double-sided hy- 
drogenation of the graphene plane. 
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